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We forecast the expected cosmological constraints from a combination of probes of both the 
universal expansion rate and matter perturbation growth, in the form of weak lensing tomography, 
galaxy tomography, supernovae, and the cosmic microwave background incorporating all cross- 
correlations between the observables for an extensive cosmological parameter set. We allow for non- 
zero curvature and parameterize our ignorance of the early universe by allowing for a non- negligible 
fraction of dark energy (DE) at high redshifts. We find that early DE density can be constrained 
to 0.2% of the critical density of the universe with Planck combined with a ground-based LSST-like 
survey, while curvature can be constrained to 0.06%. However, these additional degrees of freedom 
degrade our ability to measure late-time dark energy and the sum of neutrino masses. We find that 
the combination of cosmological probes can break degeneracies and constrain the sum of neutrino 
masses to 0.04 eV, present DE density also to 0.2% of the critical density, and the equation of state 
to 0.01 - roughly a factor of two degradation in the constraints overall compared to the case without 
allowing for early DE. The constraints for a space-based mission are similar. Even a modest 1% dark 
energy fraction of the critical density at high redshift, if not accounted for in future analyses, biases 
the cosmological parameters by up to 2a. Our analysis suggests that throwing out nonlinear scales 
(multipoles > 1000) may not result in significant degradation in future parameter measurements 
when multiple cosmological probes are combined. We find that including cross-correlations between 
the different probes can result in improved constraints by up to a factor of 2 for the sum of neutrino 
masses and early dark energy density. 



I. INTRODUCTION 

The measurement of the low-redshift expansion his- 
tory through observations of SNIa P, 0, combined 
with cosmic microwave background (CMB) observa- 
tions of flatness [1, Q and large-scale structure (LSS) 
measurements of a sub-critical matter density [5|-l7| 
provide strong arguments for a currently accelerating 
phase in the expansion of the universe. The simplest 
solution to these observations is provided by a cos- 
mological constant, or a uniformly smooth vacuum 
energy, with a pressure that is the negative of the en- 
ergy density. 

Among alternatives to the cosmological constant, 
the most popular are scalar field models with poten- 
tials tailored to give rise to late-time acceleration and 
current equation of state for the dark energy, w, close 
to -1 |8l4la|. These models are fine-tuned to have dark 
energy dominate today, just like the cosmological con- 
stant. However, the requirement w > —1 currently, 
does not imply that dark energy was sub-dominant 
at earlier times, specifically redshift z > 2, where we 
have no direct constraints. Even in single scalar field 
models, one could have a wide range of behavior for 
w{z), or equivalently the evolution of the dark en- 
ergy density with redshift. For example, oscillating 
w{z) models provide an example where the dark en- 
ergy density is not negligible in the past, while at the 
same time the fine-tuning of the potential is benign as 
compared to anthropic considerations [3] ■ 

Observationally, we have no direct constraints on 
the expansion rate of the universe at z > 2. One 
way to parameterize our ignorance is to allow for early 



dark energy. Such a parameterization also allows us 
to estimate the changes in the growth of structure 
compared to a ACDM model. 

Our aim in this paper is three-fold. 

1. We explore how well the present dark energy 
density and its equation of state may be con- 
strained using multiple probes that are sensi- 
tive to the growth of structure and expansion 
history. We include cross- correlations between 
these different probes and isolate their effect on 
parameter constraints and degeneracies. 

2. We include non- negligible dark energy density 
at early times (z > 2), ask how well this con- 
tribution can be measured, and how it affects 
the inferred values for the present dark energy 
equation of state w and density fldo. 

3. We discuss how the inclusion of early dark en- 
ergy affects the determination of the sum of the 
neutrino masses in spatially fiat and non-fiat 
cosmological models, along with the expected 
degradation from unknown reionization history. 

With respect to the second point above, we note that 
a study of expected constraints from baryon acous- 
tic oscillation (BAO) observations has shown that the 
presence of early dark energy can significantly bias 
the inferred values of late-time dark energy param- 
eters We may attempt to resolve this bias by 
using a calibration parameter, similar to that for SNe 
observations. However, this would lead to more than 
a factor of two larger degradation in BAO dark energy 
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constraints [17|, which stresses the importance of in- 



cluding this uncertainty in parameter constraints and 
on measuring dark energy contribution to the expan- 
sion rate at z > 2. 

At low rcdshift, the expansion rate in a model with 
EDE is designed to masquerade that of models with 
late-time {z < 1) dark energy [l^, thereby severely 
limiting the constraining ability of probes exclusively 
sensitive to the universal expansion, such as SN mea- 
surements. The handle on early dark energy must 
therefore come from a combination of expansion rate 
and matter power spectrum measurements. 

For a universe with dark energy to produce the mat- 
ter perturbations seen in the present universe, more 
structure needed to have formed at earlier times than 
for a universe without dark energy. Similarly, a uni- 
verse with a non- vanishing amount of dark energy at 
early times requires more structure to have formed 
at earlier times than a universe with only late-time 
dark energy. With a fixed large-scale normalization of 
the power spectrum, this implies a suppression of the 
matter power spectrum on small scales. 

Another effect on small scales is due to massive neu- 
trinos. These relativistic massive neutrinos can stream 
out of high-density regions on scales below the free- 
streaming length scale, and therefore suppress the per- 
turbations in the small-scale matter spectrum. Both 
the free-streaming of neutrinos and the suppression 
of small-scale power due to dark energy are late-time 
effects. However, their effects may be disentangled 
in either the CMB lensin g p ower spectrum [l9l . l20j , 
shear power spectrum |2ll423l | , or galaxy power spec- 
trum [2J-I23] (moreover see Ref. [25 for a combined 
lensing analysis). We revisit this issue including the 
presence of early dark energy because we expect signif- 
icant degeneracy between the effects of massive neu- 
trinos and early dark energy [29| . 

Previous efforts to constrain models of early dark 
energy have focused on measurements of the CMB 
temperature and polarization fields, large and small 
scale structure formation, Lyman-g forest, Gamma- 
Ray Bursts, and SNe [30l - [35| . In recent work, the 
influence of weak gravitational lensing observations 
of galactic sources at low redshifts and the CMB 
as a large redshift source has also been investi- 
gated [l^ [3^. In this work, we explore how well 
a generic model for dark energy that does not van- 
ish at early times can be constrained in conjunction 
with massive neutrinos from a comprehensive array of 
next-generation weak lensing, galaxy, CMB (includ- 
ing lensing), and SN observations, including all cross- 
correlations between the non-SN probes. We include 
SNe to help break the degeneracies with the late-time 
dark energy parameters. It is worth noting that the 
comprehensive inclusion of cross-correlations is new 
even for the case with just late-time dark energy. 

We take our fiducial cosmological model to be flat 
with a nonzero fraction of dark energy at high red- 
shift, rip = 0.01, as discussed in the next section. We 



assume that the dynamics is described by a canonical 
scalar field. In agreement with WMAP data [1, 
we set fldo = 0.742, wq = —1.0, (present dark en- 
ergy density and equation of state), il,ch^ = 0.11 (den- 
sity of cold dark matter), n^h^ = 0.0227 (density of 
baryons), il^h^ = 1.81 x 10"'^ (density of neutrino), 
ilk = (flatness), n,, = 0.963 (tilt of power spectrum), 
d?7,,,/dlnfc (running of tih), A|. = 2.21 x lO'^ 
(normalization of power spectrum), r = 0.087 (optical 
depth to Thomson scattering), Nes = 3.04 (effective 
number of neutrinos). 

The fiducial value of il^h^ corresponds to ^ nii, = 
94 fti^h^eV = 0.17 eV. For the inverted hierarchy, this 
value for the sum of neutrino masses actually implies 
2 mass eigenstates with masses about 0.065 eV each 
and a lighter cigenstate. We have approximated this 
as 2 massive neutrino eigenstates (with masses 0.085 
eV each) and one massless neutrino eigenstate. We 
included the fraction of ^ rrii, in one of the two mas- 
sive eigenstates as a parameter f^i but found that the 
Fisher matrix constraint on that fraction to be ~ 1 /2 
with all the cosmological probes and parameters in- 
cluded. Therefore, our results are not expected to 
change if, for example, we split ^ rtii, among 3 mass 
eigenstates. For brevity, we do not show this fraction 
parameter in the tables and figures. Even though cos- 
mology is primarily sensitive to the sum of neutrino 
mass eigenstates, it can probe the neutrino mass hi- 
erarchy in the sense that if the sum of the neutrino 
mass is below 0.095 eV (see Fi gJTT]) . then the normal 
hierarchy is the only solution |37|. 

It is important to note that our cosmological model 
above includes information about reionization only 
through the optical depth. This is sufficient if reion- 
ization happens rapidly enough compared to the age of 
the universe during reionization. However, the reion- 
ization process could be more gradual or more com- 
plicated (for example, occur in two stages) and this 
would manifest itself in a range of effects on the large- 
angle CMB polarization even when the optical depth 
is fixed [13. As a result, the estimate of how well r 
can be measured is dependent on the features of global 
reionization [38l - l40| . 

The constraint on r determines how well the am- 
plitude of the primordial power spectrum is measured 
such that (T(A|j.)/A|j = 2cr(T), where <j{X) denotes 
the uncertainty in parameter X. For constraints on 
parameters such as neutrino mass, the knowledge of 
the overall power spectrum normalization is impor- 
tant [1^. Here, we test the effect that including a 
detailed reionization model would have by imposing 
a floor on t as advocated in Refs. (T9I . [391 140| . The 
main results we quote in this paper, however, assume 
a sharp reionization scenario because our main aim 
is to see how adding different probes helps break de- 
generacies in the presence of an unknown dark energy 
component at high redshift. 

We describe our method to include early dark en- 
ergy in Section 2, reviewing the nonlinear matter 
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FIG. 1. Energy density (top) and equation of state (bot- 
tom) of early dark energy and a cosmological constant. At 
low redshifts the EDE mimics a dark energy component 
with the same density and EOS at present, and decouples 
after redshifts of a few, the exact redshift depending on 
the size of the EDE fraction Qe. 



power spectrum, weak lensing power spectrum, galaxy 
power spectrum, SN distances, and the CMB. In Sec- 
tion 3, we provide prospective early dark energy and 
sum of neutrino mass constraints (and potential bi- 
ases) obtained from a joint analysis of these probes, 
based on a Fisher matrix prescription. Section 4 con- 
cludes with a discussion of our findings. 



II. METHODOLOGY 

We begin with an overview of our calculation. We 
briefly describe the EDE cosmological model, and 
then discuss the relevant observational variables. 



A. Early Dark Energy 

Early dark energy changes the expansion rate and 
hence cosmological distances. It also changes the 
growth of density fluctuations in the universe and 
hence the matter power spectrum [T6j . 

We begin by expressing the expansion rate of the 
universe in terms of the dark energy density ^ld{z) (in 
units of the critical density) as 

, „ /a-(i + z)^ + r!„,(i + z)3 + + zf 

'^^'^ = ''"V Y^^) ' 

(1) 

where = 100 /i km s ^ Mpc is the Hubble 
constant, and {fir, ^m, ^k} are the present radiation, 
matter, and curvature densities in units of the critical 
density. The present matter density is further com- 
posed of the densities of the cold dark matter, baryons. 



and massive neutrinos {flm = ^c + ^b + ^u)- The evo- 
lution of dark energy is conventionally expressed as a 
fimction of its equation of state (EOS), w{z), 

f^rfpgg / r , i + w{z) \ 

f^.(z) = ^exp(^3y^ dz^^), (2) 

where r^^o ~ f^d(O) is the present density of the dark 
energy. 

A uniform and constant vacuum density (w ~ —1) 
is simple but suffers from the well-known coincidence 
problem. The value of the dark energy density has to 
be fine-tuned so that it only affects the dynamics of 
the universe at present. This coincidence problem mo- 
tivates the exploration of solutions other than ACDM 
(e.g. [13, EH)- Among the possibilities that allow for 
w > —1 are models in which the evolution of the dark 
energy density is such that it is large enough to affect 
the universal dynamics even at z > 2. They may even 
alleviate the coincidence problem [l^ \^ . 

A realization of early dark energy is given by 
the "tracker" parameterization of Doran & Robbers 
(2006) [l^l, where the dark energy tracks the domi- 
nant component in the universe. For this case, it is 
simpler to parameterize the dark energy density evo- 
lution directly, rather than express it in terms of an 
evolving equation of state. We use a modified form 
of the Doran and Robbers (2006) [l^ parameteriza- 
tion that tracks the equation of state of the dominant 
energy, as shown in Fig. [TJ 

n 4- 7^3+3u)o 

o,(«)^o..L_)_ 

+ nMz){i-^-j^), (3) 

hl{z) = ndoil + zf+^"'° + n^{i + zf 
+ Qrii + zf + n,,ii + z)^, 

where wq ~ w{0). The function v{z) should have the 
properties that it asymptotes to unity at large redshift 
and w(0) = 0, thus ensuring that ftd{z) asymptotes to 
fie at large redshift and f2ti(0) = Qdo- We use v{z) = 
1 — (l-|-z)'^"'° [l^, but any other parameterization such 
that dln{v) / dln{z) = 0(1) will give similar results. 
Note that the first term proportional to fldo is dark 
energy density as a function of redshift for a model 
with present density of dark energy Cldo and constant 
EOS Wo- Thus, in this parameterization with early 
dark energy, the effect at low redshift is the same as a 
model with constant EOS model wq and density Sl^o- 
Quantitatively, the fie term ("early dark energy") in 
Eqn. El constitutes [0,2.1,8.0,17.7]% at redshifts z = 
[0,1,2,3] respectively of the overall amount of dark 
energy fld{z) for w = —1 and Sle = 0.01. 

We may compute the EOS using the expression 
^(,) ^ _i + (i+^)^in[a.(3)g-(z)] _ At z=0, w{z) = ^0 
and increases with z, tending to if the dominant com- 
ponent of energy density is due to pressureless matter 
and an example of this behavior is shown in Fig. [T] 
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FIG. 2. Matter power spectrum P{k) ((Mpc//i)^) against 
wavenumber k (/i/Mpc) at z = [0, 1, 5, 10] (high to low) in 
four distinct cosmologies: ACDM (dotted black), wCDM 
(dashed red, wq = —0.9 and fie =0), ACDM with massive 
neutrinos (solid green, "^rrii, = 0.3 eV), and a CDM uni- 
verse with EDE (dot-dashed blue, fie ~ 0.05, wo = —1). 



B. Matter Power Spectrum 

Given a particular parameterization of the EOS, the 
growth history depends on the underlying mierophys- 
ieal model for dark energy. In this paper, we will 
assume that the underlying model is a single scalar 
field. Then each fid(-z) or equivalently w{z) can be 
connected to an underlying scalar field potential and 
equations for perturbations in matter, radiation and 
dark energy may be written down (e.g. [1, El].) 
This allows the calculation of the growth of matter 
perturbations. In order to correctly calculate the 
linear matter power spectrum, we numerically solve 
the Boltzmann equations with a modified version of 
CAMB [zil. We approximate the effect of an EDE 
component on the linear matter power spectrum by 
use of the PPF module by W. Fang [H, \^ . 

We extend this power spectrum to nonlinear scales 
by calculating the appropriate effective spectral in- 
dex, effective spectral curvature, and nonlinear scale, 
and employin g th e fitting functions provided in Smith 
et al. (2001) [43]. The underlying cosmology in the 
Smith et al. fitting function manifests itself in two 
distinct ways. First, cosmology impacts the evolution 
of the matter density, f2m(z), and the evolution of the 
growth of matter perturbations, D{z). The cosmology 
then fixes the functional form and coefficients associ- 
ated with the fitting functions, which are fine-tuned 
to a suite of ACDM N-body simulations. 

Thus, whereas an arbitrary dark energy EOS could 
leave its imprint on the nonlinear matter power spec- 
trum via its influence on rtm{z) and D[z), the cos- 
mological dependence of the N-body fitting functions 
remain fine-tuned to a w = —1 dark energy EOS. 
For this reason, the nonlinear solution for a non- 
ACDM cosmology is only approximate, and studies 



have shown that it could lead to an underestimation 
of dark energy constraints: for example, in prospective 
weak lensing measurements by factor of two [H, . 
Our use of this nonlinear extension could also under- 
estimate the suppression of the matter power spec- 
trum on small scales due to massive neutrinos (e.g. 
Ref [13]). Our constraints on neutrino mass, in par- 
ticular from weak lensing, may therefore be conserva- 
tive. Further work is needed to quantify this effect to 
the level of precision required for the next-generation 
telescopes considered in this paper (however, also see 
Refs. [53-i3). 

Fig. [5] illustrates the matter power spectrum in four 
distinct cosmologies: ACDM, wCDM, ACDM with 
massive neutrinos, and a CDM universe with EDE. 
The inclusion of massive neutrinos induces a sup- 
pression of the matter power spectrum on scales be- 
low the free streaming length and this is, to an ex- 
tent, degenerate with the suppression introduced by 
early dark energy. The suppression due to the pres- 
ence of early dark energy is evident in Fig. [2] and is 
a result of the fact that the universe is expanding 
faster compared to a model with no early dark en- 
ergy. The present (17do) and early {^1^ ^ 0) dark en- 
ergy densities are further degenerate with the present 
DE equation of state wq. A combination of large 
and small scale probes are therefore needed to break 
the degeneracy between these cosmological parame- 
ters {m,,^do,WQ,n,) ^MMM- 



C. Weak Lensing Tomography 

The images of distant galaxies are gravitationally 
lensed by matter inhomogeneities along the line-of- 
sight. In the weak lensing regime these percent-level 
magnifications and shape distortions of galaxies need 
to be analyzed statistically (see [s^ [s^ for a review) . 
By extracting the shear power spectrum of weakly 
lensed sources [57l - l65| . the nature of the dark energy 
has been constrained with lensing surveys [58l.[59j. An 
important aspect is that the lensing power spectrum 
depends on both the lensing kernel and the growth 
of perturbations, making lensing the most powerful 
probe of the underlying cosmology 



We employ weak lensing tomography, as this pro- 
vides information about the redshift distribution of 
the intervening lenses, and thereby allows for more 
stringent constraints on cosmological parameters 1671 
68 1 . We work within the Born approximation [63- 
7l| to compute the lensing potential (weighted pro- 
jection of the three-dimensional gravitational poten- 
tial [zHii): 



dzi^(z)Cr(z)$(n,z), 



(4) 



where K{z) = H{z)dx/dz = 1 in a flat universe. The 
lensing weight function (see Fig. of the i*^ tomo- 
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FIG. 3. Redshift dependent geometric weight Wij{z) = 

^{n,^HS)'H{z)x''^^iz)C{zKH^) (see text following 
Eqn.[T3l) of CMB weak lensing (CMBL) and tomographic 
low redshift weak lensing bins for a flat CDM model with 
= 0.1 and wo = -1.0. For CMB lensing and the fifth 
tomographic bin we also plot the kernel for Sle = 0. The 
lensing kernel captures approximately the redshift depen- 
dence of the integrand for the lensing power spectrum 
(including that from the matter power spectrum). The 
diminishment of the kernels for f2e > stems from the 
increase of H{z) with increasing Q,e- 

graphic bin is given by 

where X{zo, z^) is the comoving distance to object at 
redshift Zg measured by observer at redshift Zq, such 
that x{z) = -'^(O, 2), and Sk{x) = [sin(a;), sinh(a;), x] 
for a [closed, open, flat] universe, respectively. We set 
Ci{z) = explicitly if z > z^+i. The quantity p{zs) 
contains the source galaxy distribution, the integral 
over which is fii (see Sec. lIII A|) . For the lensing of the 
CMB photons, we replace the galaxy source distribu- 
tion with a (5-function at the last scattering surface. 

The cross correlation of the respective Fourier coef- 
ficients at angular multipoles 1 and 1' is given by 

{cb*mi^')) = i2nrSu{l-r)Ct/i£) (7) 

in the flat sky limit, which through the Limber ap- 
proximation reduces to [73l - [75| : 

Cffie) = dzK{z)C{z)Q^{z)^P^^ (fc, z) . 

Jo X [z) 

(8) 

This projected spectrum remains the same in the all 
sky formulation [72|. 

To find the relation between the power spectrum 
of the potential (P$$(k, z)) and matter perturbations 
(PssO^jz)) on sub-horizon scales, we use the Poisson 



equation in Fourier space, 

a>(k, z) = -^n^ (^^) 5(k, z){i + z). (9) 

Accounting for the definition of two-point correla- 
tions of the potential and density fields: 

($*(k, z)$(k', z)) - (27r)3fc(k - k')P*$(A:, z), 
(r(k,z)(5(k',z)) = (27r)3fc(k-k')PA-A-(fc,zXlO) 

we obtain the power spectrum of the potential in 
terms of the power spectrum of density fluctuations 
as 

P^^{k,z)^^-^ll^{l + zf{^^ Pss{k,z). (11) 

Reexpressing the power spectrum of the matter as a 
dimensionless quantity, 

Ajs{k,z) = k^Pss{k,z)/{2TT^), (12) 

the spectrum of the convergence field, Cf* = 
{1/A)e{e + ifCf'l', is given by: 

C5« (£) = AT dzK{z)C:{z)q{z) X (13) 

Jo 

(l + z)2i/(z)x^(z)Al,(^±i/?,z), 

where the density and scale dependent amplitude 

= (9^V8)(^^mHo')'^'(^+l)'(^+l/2) and 
A^^(£/x(z), z) encapsulates the full nonlinear dark 
matter power spectrum. We construct a geometric 
weight Wij{z) that approximately captures the con- 
tribution to the convergence power spectrum from 
different redshifts for a broad range of multipoles 
by removing a factor of (1 + z^x'I'^ij.') from the 
integrand in Eqn [T31 This results in Wij{z) = 
^{n,„H^)^H{z)x'^^^{z)Q{z)q{z). The reason this 
approximation takes into account the redshift contri- 
bution to the convergence power spectrum is that the 
matter power spectrum A^^(fc,z) scales roughly like 
^3/2 ^Qj. gcales ©(O.l/Mpc). We plot this geometric 
weight Wij{z) in Fig. [3l and this shows that most of 
the contribution to low redshift weak lensing comes 
from z w 0.5, and in the case of CMB weak lensing 
z « 1. 

Fig. |4] further illustrates the effect of varying the 
dark energy and other parameters on the correlation 
function in the third tomographic bin. The shapes are 
roughly the same across the other bins. The similarity 
in the shapes of the variation due to early dark energy 
and sum of neutrino masses implies that there will be a 
degradation in the constraints on the sum of neutrino 
masses due to the addition of early dark energy. 



6 




10 



1000 



10 



100 



100 

I i 

FIG. 4. Left: Top: Convergence power spectra + l)Cij"' /2ti for the case of five tomographic bins in the fiducial cos- 
mology. We include the expected noise for LSST as a band about the curve. Mid: Logarithmic derivative dlnCij'^ / d\npu 
of the convergence power spectrum with cosmological parameters for tlie third tomograpliic bin (i = j = 3). The 
derivatives of the other tomographic bins have similar shapes. Bottom : The sub-window zooms in on the logarithmic 
derivatives with sum of neutrino masses (dotted bhie) and EDE density (dot-dot-dashed red). Right: Same as Left but 
for galaxy power spectra. The noise contribution (both LSST and JDEM) is at most on sub-percent level and decreases 
towards smaller scales, as the noise is constant with scale whereas increases with scale. This is to be contrasted with 
lensing tomography where CI"" peaks at around I = IQ and thereafter decreases continuously. 



D. Galaxy Angular Power Spectrum field i5^(fi, z) is 



The galaxy-galaxy angular power spectrum C^^{£) 
can be obtained analogously to the derivation of 
C'^'^{£) in Eqn. [H] The two-dimensional projection 
of the three-dimensional biased tracer of the density 



5.(fi)= / "'^^ dzK{z)C!iz)6Sin,z). (14) 
Jo 

Froii-i{g*(k,z)g{k',z)) = (27r)3,5D(k-k')P<,g(fc, z), we 
thereby obtain the projected galaxy power spectrum 
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FIG. 5. Left: Top: Tomographic galaxy-Iensing correlations l{l -\- l)Cl^ /2ti. Here we only illustrate a subset of cases 
where galaxy and lensing bins fully overlap. Mid: Logarithmic derivatives of the power spectra dlnC^? / dlnph with 
parameters pu for the third tomographic bin in both galaxy and lensing (similar characteristics for other bins). Bottom : 
Zooming in on the derivatives with EDE density and sum of neutrino masses. Right: Logarithmic derivatives of the CMB 
temperature (solid black), E-mode (dashed blue), and lensing potential (dotted red) power spectra with cosmology. 



in the Limber approximation: 

Jo X \^) 

(15) 

where the weights for Zi < z < Zi+i are 

Cf(z) = p(z)M, (16) 

and C,f{z) = if z ^ (zi,Zi_|_i). As a result, only the 
auto-correlations survive the integral. In the linear 
regime, the three-dimensional galaxy power spectrum. 



Pgg{k, z), is related to the matter power spectrum via 
the bias b{k, z): 

Cffii) = Af dzK{z)C!{z)Cl{z) (17) 

Jo 

X H(.)x(z)62(A:,z)A2, (^^,^-), 

where amplitude = 27r^ {i+ 1/2) ^ . Restricting 
the analysis to largely linear scales (see Table |ll|) , we 
include a time-varying linear bias by the first-order 
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FIG. 6. Left: Top: Correlations between weak lensing of the CMB and low-redshift sources ^(^ + l)Cf^''/27r for three of 
the five tomographic bins. Mid: Logarithmic derivatives of the power spectrum din C,''^'" /din with parameters for 
the third tomographic bin (similar characteristics for other bins). Bottom : Zoom on derivatives with EDE density and 
sum of neutrino masses. Right: Same as Left but for correlations between galaxy tomography and lensing of the CMB. 



expansion b{z) = 5o + biz. A Gaussian prior of 10% 
is applied to these two bias parameters [76l 1771 . and 
the fiducial values are set to [fco,^i] = [1.0,0.8] [T^I. 

E. Galaxy-Lensing Correlation 

The cross-correlation of the source galaxy distribu- 
tion and the weak lensing convergence is: 

Jo X [z) 

(18) 



Correlations between tomographic bins where the 
galaxy bin lies behind the lensing bin are zero (result- 
ing in 15 out of 25 non-zero correlations). The power 
spectrum of galaxy and potential perturbations may 
be written as: 

Pg$(fc,z) = '^^^(jY^ b{k,z)rc{k,z)Pss{k,z), 

(19) 

where we have introduced the correlation coeffi- 
cient Tc defined through the relation Pgg{k,z) = 
rl{k,z)Pgg{k,z)Ps5{k,z). Then, using Cf = + 
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FIG. 7. Left: Top: ISW-WL correlations + l)Cf "^/27r for three of the five tomographic bins. Mid: Logarithmic 
derivatives of the power spectrum dlnCf "^/dlnpfc with parameters for the third bin (similar characteristics for other 
bins). Bottom : Zoom on the EDE and neutrino mass derivatives. Right: Same as Left but for ISW-galaxy correlations. 



1)C^ /2, we may write the angular power spectrum 
for galaxy and weak lensing correlation as: 

Cf^ii) = Ar dzK{z)C!iz)giz) X (20) 

JO 

(1 + z)H{z)x'{z)b{k, z)r,{k, z)Ajs (^^' ^) ' 

where we have defined = -{Sn^ /2){nmH^)i{i + 
l)(i' + l/2)~^. Restricting the analysis to linear scales, 
we fix = 1 in our calculations. 



F. Cosmic Microwave Background (T, E, Kc) 

We consider CMB information from the tempera- 
ture field (T), E-mode polarization (E), and weak 
lensing of the CMB (kc)- The temperature and E- 
mode polarization power spectra (and their cross- 
correlation) are obtained from a modified version of 
the Bohzmann code CAMB 

As for lensing of the CMB photons [Hi 113, [HI (also 
see Refs. [H, |83 for detections of the lensing power 
spectrum), we calculate it in the same manner as in 
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Egn.lTSlwith the source constrained to be at the red- 
shift of decouphng (zdec), such that 



H{z)x{z) Xizdec) 



(21) 



Hence, aside from Q{z)C,'^{z) — !• C,'^''{z)C,'^''{z) (and 
^max = 2^dec), C"*'='''=(£) has the Same functional form 
as Cij^(£). Similarly, the cross-correlation spectrum 
C^^-'ll) has the same form as C^'-'{£) with the sub- 
stitution Cj{z) C''iz), and we obtain '=(^) by 
replacing Cj{z) with C'^^{z) in Cf^[£). 

In Figure [5] we show the logarithmic derivatives of 
the CMB temperature, E-mode, and lensing potential 
power spectra with cosmological parameters. In par- 
ticular, it is evident that both CMB temperature and 
lensing measurements of the early and late universe 
are necessary to break the degeneracies between early 
dark energy and sum of neutrino masses. 



G. ISW-Lensing Correlation 

Wc consider five distinct lensing cross-correlations: 
the first is that of weak lensing in redshift bins {i) 
with the lensing of the CMB photons, Cf'''=(£); the 
second is that of galaxy counts in redshift bins {i) 
with the lensing of the CMB photons, Cf'^''{£); the 
third is that of galaxy counts in redshift bins (i) with 
weak lensing in redshift bins (j), Cf^^{i); the fourth 
is that of weak lensing in redshift bins (i) with the 
unlensed temperature, Ci'^{i); and the fifth is that 
of lensing of CMB photons with the unlensed tem- 
perature, C'^"'^ {€). Having already described the first 
three correlations, we here discuss the last two corre- 
lations that arise due to the Integrated Sachs Wolfe 
(ISW) effect ll-il. 

The perturbations in the CMB temperature due to 
the late-time gravitational redshifting of the photons 
is encapsulated in [t^] 



T 



dzK{z)C{z)Mi^,z), (22) 



where the overdot denotes a derivative with time t, 
and the ISW weight is 



{l + z)H{z) 



(23) 



The cross spectrum of the respective Fourier coeffi- 
cients is given by 

i(0:(l)Ar(l')) = (2^)2fc(l-l')CrW, (24) 
which via the Limber approximation reduces to 



Cf{l) = / dzK{z)Q{z)e{z)^P^^{k,z). 
Jo X (Z) 

(25) 



To find P^^{k, z), the time derivative of the potential 



2(l + z)-i \ k J \H-^{z) 



and hence 



(26) 



^$$(^1 z) ^ . 



9 nl 



4(l + z)-2 y k 



Pss 



H-\z) 



(27) 

To compute Pj^, we note that in the linear regime 
Pgg = ^^Pss- This relation has been numerically 
verified to hold in the nonlinear regime up to £ = 5000 
[sst . Hence, we adopt this relation so that 



^tjxi) (^j z) 



9 nl^Hjz) (HoVf^ , 1 



4(l + z)- 



+ z d 



2 dz 



Pss, 
(28) 

and the resulting k — T angular power spectrum is 
then given by 

Cr(£) = Af r^\zK{z)Q^{z)C{z)x'{z) x (29) 
Jo 

1/2 



(l + z)^i/^(z)(l + i±i|)A^, 



X{z) 



where the amplitude Af (£(£+l)/2) ^ A'l'^ = 
(97rV4) {nmH^)^e{£ + 1){£ + 1/2)-^. As the deriva- 
tive I'^^^^aal < and the lensing weight con- 
tains a negative sign, it follows that C^'^{£) < 0. We 
obtain C'^'''^{t) by replacing C,i{z) with C,'^''{z) in the 
above equation. 



H. ISW-Galaxy Correlation 

To obtain the ISW-galaxy correlation, 

E{z 
'X^Y'^' 



Cf{£) = I dzK{z)(:f{z)e{z)^^P^^{k,z), 

(30) 

we follow the same steps as for the ISW-lensing cor- 
relation and obtain 

P,i{k,z)^\^^i^^^\{k^z)r^{k,z) (31) 
X {H{z)Pss{k,z)-P,^{k,z)). 
We thereby find the projected spectrum 

Cf{£) = Af dzK{z)C!{z)e{z) (32) 

Jo 

X {1 + z)H^{z)x'{z)b{k,z)T,{k,z) 



X 1 



1 + z a\ 2 ^^ + 1/2 



,z \ , 



2 dz) \ x{z) 
where amplitude Af = -37r2(17,„iJ^)(f + 1/2)"^ 
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I. Supernovae 

We include the supernovae distances in our analysis 
for their ability to constrain the expansion history of 
the low-redshift universe [89j . The observable quan- 
tity is the apparent magnitude: 



5 log 



10 



(i + ^)x(^) 

Mpc 



25 + M, 



(33) 



where M is the absolute magnitude, and we marginal- 
ize over it around the fiducial value of Mfid = —19.3 
with a Gaussian prior of 0.6. This is equivalent to 
treating the distance modulus ^{z) = m{z) — M as the 
observable, and adding to that the nuisance parame- 
ter M. . The fisher matrix of the prospective supernova 
measurements is then 



dzN{z) 



dm dm 

dpa dpj3 ' 



(34) 



where am is the measurement uncertainty, N{z) is the 
redshift distribution of the SNe, Zmax is the redshift 
of the farthest SN, and p represents the cosmological 
parameters. 

Although, we do not expect a strong measurement 
of EDE from SNe, we will see that it does help in 
breaking some of the degeneracies in the parameters 
of the background cosmology, in particular the dark 
energy at the present epoch. To this end, we uniformly 
distribute a set of 300 SNe at z < 0.1 [H [9|. For 
the space-based JDEM survey, we add 2000 SNe in 
the range 0.1 < z < 1.7 [U, whereas for the LSST 
survey we distribute 3.0 x 10^ SNe (corresponding to 
six years of data) between 0.1 < z < 0.8 For each 
supernova we take the intrinsic noise to be Gaussian in 
magnitude with (Tint =0.1 [9^. We divide the Hubble 
diagram for z > 0.1 into 50 redshift bins, and associate 
each bin with a redshift-dependent systematic error 
floor of magnitude 

Smiz) = 0.01(0.1/Az)i/2(i.7/^„^^^)(l + 2)/2.7, (35) 

where Az is the width of the relevant redshift bin [9^ . 
For the ground-based LSST survey, we allow for pho- 
tometric uncertainties of the form 



(Tpjl + z) 



(36) 



where iVc = 100 is the number of spectra used in the 
photo-z calibration in each bin, and ap = 0.01. As- 
suming no correlation between bins, the total magni- 
tude uncertainty is thereby quantified via 



<ymiz) 



1 



fdm\ ^ 



V d 



Sliz), (37) 



where iVbin is the number of SNe in each redshift bin, 
and (Tj = for the spectroscopic JDEM survey. The 



dominant contribution to the total error comes from 
the systematic floor, and the photometric redshift er- 
rors turn significant only for more pessimistic choices 
of ap- 

For a fixed Zmax, the SN constraints on the under- 
lying cosmology are robust to variations in the par- 
ticular SN redshift distribution of the survey, while 
this maximum redshift itself imposes an important 
limiting factor on the cosmological constraints [9^ 
[97| . Due to the systematic fioor, we find the LSST 
constraints on the underlying cosmology from SNe re- 
main effectively unchanged after the first year's data. 
The existence of a systematic floor combined with SNe 
extending to higher redshifts explain why our JDEM 
configuration provides stronger constraints than our 
LSST configuration when it comes to SN explorations 
of dark energy (as shown in Table IIVI) . 



III. RESULTS 

We next explore the constraints on early dark en- 
ergy from future weak lensing, large-scale structure, 
SN, and CMB measurements (and all relevant cross- 
correlations). We begin by listing the observational 
properties of these probes. Then we elaborate on the 
constraints on the underlying cosmology of the uni- 
verse, and the biases in the cosmological parameters 
that arise due to a potential neglect of early dark en- 
ergy. 



A. Survey Properties 

We consider Planck HI] for CMB temperature, 
polarization, and lensing measurem ents, and compare 
this to a future mid-cost CMBPol [lOO, [l03 mission 
(Epic-2m). For the large-scale structure, weak lensing, 
and SN observations, we compare two generic types of 
surveys: 1) a "wide" survey, and 2) a "deep" survey. 

For the "wide" configuration, we consider a survey 
covering half of the sky, such as the gro und-based 
Large Synoptic Survey Telescope (LSST) [loi, [l0|. 
For comparison, for the " deep" configuration, we con- 
sider a survey covering a tenth of the sky, such as 
a space-based Joint Dark Energy Mission (JDEM) 
candidate li ke the S uperNova Acceleration Probe 
(SNAP) [H [Tog - [l06t . For weak lensing and galaxy 
clustering measurements, beyond the sky coverage, we 
only allow for the number density of source galaxies to 
be different between these two survey configurations 
(as seen in Table For the SNe, the two survey con- 
figurations are different in the number and redshift 
extent of observations (as seen in Sec. Ill ip . We re- 
fer to the wide survey as LSST, and the deep survey 
as JDEM. Other surveys w ith similar c haracteri stics 
include WFIRST [io3,IiQ3 and Euclid [T09l . fTToj . 

The distribution of source galaxies is divided into 
five tomographic redshift bins, as the gain in cosmo- 
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Large-Scale Structure Survey Properties 



Probe 


/sky "-g 


(arcmin ^) 


■^pcak 




^max 


No. bins 


LSST 


0.5 


50 


1.0 


0.22 


3000 


5 


JDEM 


0.1 


100 


1.0 


0.22 


3000 


5 


LSSTiooo 


0.5 


50 


1.0 


0.22 


1000 


5 


JDEMiooo 


0.1 


100 


1.0 


0.22 


1000 


5 



TABLE I. Description of the ground-based (LSST) and 
space-based (JDEM) probes. The redshift bins of the 
source distribution are delimited at [0.75, 1.1, 1.45, 1.95, 
3.0], such that each tomographic bin contains roughly the 
same number density of galaxies. The minimum angular 
multipole ^min = 2 (see text for more details). 



logical information diminishes rapidly for more than 
five bins (49l . \q1\ . The number density of sources in 
a square arcminute in each redshift bin (with bound- 
aries Zi < Zs < Zi+i) is defined by 



Zi + l 



dzs p(zs 



(38) 



where the redshift distribution of source galaxies is 
taken to be of the form [571 ): 



(39) 



a = 2, P = l} [log, appro- 
llll| . normalized such that 



We adopt {zq = 0.5 
priate for LSST [11 

/g dz p{z) = fig. We use the same distributionto 
describe our JDEM source population [ol [lOl [l05| 
with a factor of 2 higher fig value (see Table [l|). 

The observed convergence power spectrum, identi- 
cal to that of the shear [H^, is contaminated by shot 
noise due to the finite source density, as well as un- 
certainty in the intrinsic shapes of the source galaxies, 
leading to 



^(c^rW + %(7')M). (40) 



which assumes that the noise is uncorrelated between 
tomographic bins. Similarly, for the galaxy angular 
power spectrum. 



(41) 



We take the intrinsic shape uncertainty of the source 

1/2 

galaxies to be redshift independent: (7^) = 0.22, 
in accordance with expected results for the future 
ground-based probe. For simplicity, we keep the same 
source distribution and intrinsic shear uncertainty for 
the surveys on ground and space, modifying only the 
source density to twice that of the ground-based sur- 
vey, and the width of the survey to a tenth of the 
sky. Table [J summarizes the characteristics of the two 
surveys. 



=j=2 LSST 
=j=2 JDEM 
=j=5 LSST 
=j=5 JDEM 
=2, j=5 LSST 
=2, j=5 JDEM 




FIG. 8. Weak lensing signal-to- noise Ci/'ie)/ ACtJ'il) for 
combinations of the second and fifth tomographic bins in 
the fiducial cosmology, for LSST (thick) and JDEM (thin). 



In Figure [H we plot the power spectrum of the con- 
vergence for a subset of the five tomographic bins 
(with the redshift divisions listed in the caption of 
Table [T| . We divide the power spectrum by its noise, 

^c^;ii) = y (c-r wc-j^ w + At (42) 

in Figure [5| Whereas cosmic variance dominates the 
error on large angular scales, the shot noise is dom- 
inant on small scales. The signal to noise is consis- 
tently higher for the wider LSST than it is for the 
deeper JDEM. The larger width of LSST gives strong 
signal to noise in particular at \ow-i, while the greater 
depth of a JDEM-type survey makes it increasingly 
competitive at high-^. 

For the weak lensing surveys we analyze multipoles 
between £ = 2 — 3000 in one case, and multipoles 
of ^ = 2 — 1000 as a second case. The cutoffs at 
i = [1000, 3000] la rgely avoi d non-Gaussianities of the 
convergence field [I12l4ll4 . as well as uncertainties 
from bary onic physics that increase at larger multi- 
poles jll5^1Tl8i . We impose a cutoff in the matter 
power spectrum for k > fcmax = 10 /i/Mpc, which 
implies that for the smallest angular scales, £max = 
[1000,3000], distances below x = [100,300] Mpc/h 
(or equivalently, redshifts below z ~ [1/30, 1/10]) are 
wiped out in the calculation of the convergence spec- 
trum. Thus, the cutoffs ensure that scales smaller 
than k — 10 h/Mpc in the matter power spectrum are 
only probed at low redshift, where the number densi- 
ties of sources approach zero. Recent hydrodynamic 
simulations that include AGN feedback to solve the 
overcooling problem have shown these cutoffs may not 
be conservative enough unless the baryoni c effects i n 
the power spectrum arc precisely modeled [ll7l IllSj ] . 

At the other end of the spectrum, we keep all mul- 



tipoles down to £„ 



2. We incorporate perturba- 



tions in the dark energy fluid on these scales. While 
a JDEM-type survey is not expected to probe scales 
below £ ^ 10, we have checked that our extension 
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Nonlinear Cutoffs in Galaxy Surveys 



CMB Survey Properties 



Bin 12 3 

^median 0.38 0.93 1.3 



4 5 
1.7 2.5 

fcmax 0.11 0.14 0.17 0.21 0.29 
120 320 490 720 1200 

TABLE II. We effectively eliminate nonlinear scales 
in the constraints from galaxy surveys by setting 
[Cl^,Cf, Cf/, Cf"'=] {t > £m.^) = 0. The minimum an- 
gular scale £max is computed via £ + 1/2 = kx, where we 
let X = X (^median) and approximate, for the considered 

redshifts, k = fcn,ax (Zmcdian) 0.1 X (1 + ^Cdian/S) Ell. 



down to £ = 2 does not significantly affect ou r re- 
sults. We use the Limber approximation [tI, Ill9l | 
for all angular scales for computational reasons. This 
will lead to some errors in the constraints but this is 
not a significant source of concern because the Lim- 
ber approximation has been shown to work better 
than 3% for ( > 20, better than 5% for £ > 10, 
and bet ter than 30% even for £ > 2 (see Fig. 19 
of Ref. [l20j ). For the galaxy surveys we eliminate 
nearly all nonlinear scales in our analysis by set- 
ting [G,f ,Cf ,Cf;,Cf1(^>4.ax) = 0. The min- 
imum angular scale £max is computed via £ ~ kx, 
where we let x = X (^^modian) and approximate, for the 
redshifts under consideration, k = fcmax (^^median) — 

0.1 X (l + ^i^Cdian/^) [HI- The cutoff scales are listed 
in Table HH 

In the next section (Sec. IIIIBp we will find, as also 
reported by the Dark Energy Task Force [6y|, weak 
lensing to be strongest future probe of the underly- 
ing cosmology of the universe, owing to its sensitivity 
to both structure formation and the universal expan- 
sion. However, the optimism associated with lensing 
is predicated on overcoming the vast systematic un- 
certainties in both measurement and i n theory (isl - 
Isol. [67l. [69l-[7ll. flT l flTsl fil l [T2l - [l32j . These sys- 
tematics include dark energy corrections to the mod- 
elin g of the nonlinear matter power spectrum 
[sol . Il3l| . higher order correction terms in the lens- 
ing integral (such as due to the Born approxima- 
tion and lens-lens coupling [69l - [7l| ). and uncertain- 
ties of the matter power spectrum on nonlinear scales 
due to the strong influence of baryonic physics [llSl - 
IllSl . Il29l . Il30l Il32| . Observational systemati cs include 
intrinsic galaxy alignments (e.g. (l23l4l28 [). photo- 
metric redshift uncertaint ies, shear calibration errors, 
and PSF anisotropies [67l. [l22l . [T26| . 

Furthermore, the observed ellipticities of weakly 
lensed galaxies are sensitive to the reduced shear, 
g = 7/(1 — k), where 7 is the shear and k is the 
convergence. In the weak lensing regime we make use 
of the expansion of the reduced shear to first order in 
the fields: g ~ j. For future lensing surveys, it has 
been shown that this approximation induces a bias on 
the cosmological parameters at the sam e order as that 
of the parameter constraints [l33l Il34| . Our main mo- 
tivation is to elucidate the degradation in constraints 



Experiment Channel FWHM AT/T X 10^ AP/T X lO^^ 
Planck 100 10 25 40 
143 7.1 16 30 
217 5.0 24 49 
EPIC-2m 100 8.0 0.84 1.19 
150 5.0 0.81 1.15 
220 3.5 L24 1.75 

TABLE III. Experimental specifications for the Planck 
and mid-cost CMBPol (EPIC-2m) missions. The sky frac- 
tion /sky ~ 0.65, and the angular multipoles extend from 
fmin = 2 to i'max ~ 2000. The channel frequencies are 
given in GHz, and the angular resolutions in arcminutes. 



due to EDE and the biases in the cosmological pa- 
rameters from neglect of EDE. Wc therefore continue 
with the above mentioned assumption of the shear as 
the lensing observable. 

Intrinsic galaxy alignments reflect the reality that 
underlying galaxies that become lensed are not in- 
herently circularly shap ed, but p ossess a non-trivial 
intrinsic ellipticity (e.g. |l23l - ll28| ). The assumption 
that these intrinsic ellipticities are uncorrelated with 
each other and the gravitational shear breaks down 
for high-precision measurements of weak lensing from 
next-generation experiments. For instance, nearby 
galaxies that lie in the same gravitational potential 
of a massive dark matter halo will experience the 
same tidal for ces, giving rise to an intrinsic-intrinsic 
T23H128I . ver, intrinsic-shear correlations 



term 

occur for pairs of galaxies where one galaxy is corre- 
lated with the surrounding density field in the fore- 
ground, thereby contributi ng to the lensing distortion 
of the background galaxy 123l4l28 |. 

While intrinsic ellipticity correlations can be con- 
trolled by downwcighting close galaxy pairs from the 
analysis |135. J32] , the removal of shea r -ellipticity 
correlations is more difficult [l23l - [T2l IT38Ml40l |. 
Ref. [l24| has shown that the dark energy EOS may 
be biased by 50% if intrinsic alignments are not ac- 
counted for in next-generation weak lensing exper- 
iments. The presence of intrinsic alignments also 
tighten the requirements on photometric redshift un- 
certainties by a factor of three, rega rdless of the frac- 



tion of catastrophic outliers |124| . Suggested ap- 
proaches to a void these degradations include halo 
modeling 128 1 and nulling techniques jl39l |. Alter- 
natively, the systematic may be self-calibrated for an 
extended set of observables that can correlate with 
the weak lensing shear and intrinsic alignment sig- 
nals, such as g alaxy clus terin g and lensing magnifica- 
tion dH [l40[ (also see [HI). 

The effects of multiplicative and additive uncer- 
tainties in the convergence (e.g. from shear miscal- 
ibration and PSF anisotropies, respectively) on the 
measurement of dark energy are influenc ed b y the 
correspondin g pri ors. As shown in Ref. |l22| (also 
see Refs. [i^ 114 |). the cosmological parameter con- 
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straints are compromised for multiplicative systemat- 
ics at the 1% level, and mean additive shear system- 
atics at the 10 ~^ level. The situation is analogous for 
the uncertainty in the photometric redshift distribu- 
tion of the sources, where the parameter constraints 
from Icnsing arc either heavily influenced (> 1% prior) 
or minimally influenced f< 0.1% prior) by the photo- 
metric uncertainties [i^, [63, Il22| . Fortunately, it has 
been shown that a complementary spectroscopic sam- 
ple of 10^ — 10^ galaxies efficiently protects against 
photom etric redshift errors as well as catastrophic 
outliers |l4l| . whereas alternative methods may even 
satisf y the syst ematic requirements from photometry 
alone pl[l43t . 

Thus, in this work, we will assume that these sys- 
tematic difficulties have been largely overcome with 
minimal influence on the constraints by the time 
the data from the considered next-generation lens- 
ing probes are analyzed. At the same time, we are 
not incorporating further statistics that can be ex- 
tracted from weak lensing, such as that included in 
the bispectrum [144| - |147 1 , or utilizing the complemen- 
tarit y between measurements of shear and magnifica- 
tion [l25lfl49t . 

We end this section with a summary of the CMB 
temperature, polarization and lensing noise proper- 
ties. The effective experimental noise power spec- 
trum associated with the temperature and polariza- 
tion fields is given by a summation over the number 
of channels, 



7V°°(^) 



E 

. 1=1 



Aa 



„i(i+l)6li/161n2 



(43) 



where Aa is the detector noise for a G {T,E), 6 de- 
notes the beam FWHM, and wc assume N^^{t) — 0. 
The optimal noise power spectrum of a qua drat i c es- 
timator of the convergence field is given by |l50l . Il5l| 



TT T7 



X + l)/2f {21 



(44) 



where T denotes the lensed temperature, and 



(2Zi + l)(2/ + l)(2 /2 + l) h I h 

I 



4-77 



1, 



X -[I{l + l) + l2{l2 + l)-h{h + l)], (45) 

where the quantity in brackets is the Wigner-3j sym- 
bol. Finally, we define 



^2e'+T iC'"'{i) + SabN^''{i)), (46) 



where {a, 6} G {T,E,Kc}. Values for the considered 
CMB experiments are given in Table IIIIl Secondary 



non-Gaussianities in the covariance from the trispec- 
trum (due to weak lensing, the ISW effect, and the 
SZ effect) have been shown to degrade the Planck 
and EPI C parameter constraints by 20% and 30% 
jl52lll54j respectively; however, their full account lies 
beyond the scope of this work. 



B. Comprehensive Parameter Forecasts 

In previous sections we explored the qualitative in- 
fluence of EDE on the lensing, galaxy, supernova, and 
CMB obscrvables, via its impact on the expansion rate 
and matter power spectrum. We now examine how 
these corrections quantitatively affect the combined 
constraints of the dark energy. To t his end, we utilize 
a Fisher matrix formalism 1751. Il55l : 



= ^ X Tr 



F 



■SN 



(47) 

where the decoupled SN flsher matrix is deflncd in 
Eqn. I34[ and for the combined observational analysis 
the symmetric matrix 
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(48) 

such that {k} consists of the spectra from five tomo- 
graphic bins (ki, K2, ^3, K4, K5) and {g} consists of the 
spectra from five tomographic bins (51,52753,54,95)- 
Cj''^^''\ Cj^^^^^ Cj''"^^ are therefore 5x5 subma- 

trices, and cj'"^''', cj"^"^, cj^^''', c\^^'^ are 5 x 1 
submatriccs. For the terms in Eqn. 27] we carry out 
two-sided numerical derivatives with steps of 2% in 
most parameter values. Wc have confirmed the ro- 
bustness of our results to other choices of step size. 

In Tables irvTIXl we illustrate prospective constraints 
from Planck/EPIC CMB temperature (T), E-mode 
polarization {E), lensing (kc), LSST/JDEM weak 
lensing tomography (k), galaxy tomography (5), SNe 
(s), and their combined impact (including all relevant 
cross-correlations shown in Eqn. 05]) on the 12 consid- 
ered cosmological parameters (fido, ^^e, i^ch^, ^bfi^, 
^k, NcS, Wo, n-s, dn^/dlnfc, A|, r). 

The contents of our tables are as follows: In Ta- 
ble IIVI and Table [V] we consider only a fiat universe, 
with curvature always considered in the other tables. 
These tables present the separate constraints on the 
underlying cosmology obtained from the CMB, lens- 
ing tomography, galaxy tomography, and SNe, along 
with the synergies attained from a combined analy- 
sis of these probes. Table [V] differs from Table IIVI in 
that it fixes the early dark energy density. Table I VII 
differs from Table IIVI in that it allows for variation 



15 




> 

0) 



0.000 



0.054 




0.000 



0.054 



1.05 



d" 0.96 - '*'':S^?_-> 



0.87 



2.64 
0.00 



0.000 



0.027 



0.054 



-2.64 

T3 



= , 



0.000 



0.027 



0.054 



5.4 



/ 2.7 



0.0 




0.000 0.027 



0.054 



3.06 
2.21 
1.36 







0.000 0.027 



0.054 



0.160 



0.110 

a 



0.060 
0.000 




I 

O 



411 



0.054 



3.47 
2.27 
1.07 



.MM"".,** 

. . • * 



0.000 



0.027 



0.054 



0.742 




I 

O 



0.705 



9.19 



8.70 : 



8.21 



0.000 



0.054 



/ / 


\ \ 
1 V 


- / + 


/ ; 




/ /. 


■ ''v \ 


^ 


'^1, 





0.002 



0.010 



0.018 



FIG. 9. Parameter degeneracies with early dark energy density (fie) in a flat universe (also see Table lrV|l for Planck 
measurements of temperature and CMB lensing spectra [TT, EE, TE, k.cKc, KcT] (dot-dashed, black), along with 
LSST tomographic weak lensing spectra [kk] (dotted, blue) and tomographic galaxy spectra [gg] (dashed, turquoise). 
Constraints from SNe are too weak to be visible in the shown parameter regions. The error ellipses from the combination 
of all these probes (including SNe), incorporating all cross-correlations (see Egn. I48p is shown as (solid, red) curves. 
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in curvature. In Table [Xj we present results where 
the CMB constraints are derived f rom a fut ure exper- 
iment like the proposed 2m EPIC [lOOl . Il0l| (compare 
to Table |Vl|. Table EIIl differs from Table Ell in its 
neglect of SNe measurements. Lastly, Table HXl differs 
from Tabic IVII in that we neglect cross-correlations 
between the obscrvables (i.e. neglecting all correla- 
tions between T, E, Kc, <7, except for tomographic 
cross-correlations within k). 

We now explore each of these tables in greater de- 
tail. Table HVl shows us that the dominant constraint 
on the fraction of dark energy at early times is drawn 
from the CMB (in particular TT and to some extent 
Tkc) due to its deep redshift information. At a value 
of cr(17e) = 8.6 X 10-3 (Table HyI, the Planck CMB 
temperature and polarization constraint is within a 
percent of the critical density. In general, EDE is 
best constrained by the CMB, followed by weak lous- 
ing tomography, galaxy tomography and SNe in that 
order. For comparison, the low-redshift lensing con- 
straint from LSST on fig is a factor of four (factor of 
six for JDEM) weaker than from the CMB. If we im- 
pose a nonlinear cutoff to the convergence spectra at 
^max = 1000, the situation becomes more dire, as the 
LSST and JDEM lensing constraints become worse by 
another factor of three. Similarly, the galaxy tomog- 
raphy constraint from LSST on fig is a factor of five 
(factor of nine for JDEM) weaker than the CMB con- 
straint, and the LSST SN constraint is a factor of 70 
weaker (factor of 30 for JDEM) than the CMB. 

Nevertheless, once the six observables (T, E, Kc, k, 
g, s) and all relevant cross-correlations (see Eqns. l47l 
SHI) from Planck (or EPIC) and LSST (or JDEM) are 
analyzed in a combined setting, the constraint on fie 
improves by a factor of four over the CMB constraint. 
The combined constraints are equally strong regard- 
less of the choice of LSST or JDEM for the non-CMB 
observations (k, g, s). For a JDEM-like experiment, 
the cross-correlations improve the fig constraints by a 
factor of about 2. 

As expected, we find the late-redshift parameters 
more strongly constrained by the non-CMB probes. 
For example, in a universe where we allow for the ex- 
istence of early dark energy, the LSST weak lensing 
constraints on the present DE density (ildo) and EOS 
{wq) of 1% and 6% are much better than the con- 
straints obtained from just CMB lensed data of about 
7% and 20% on present DE density and EOS (Table 
lIVp . Galaxy tomography measurements with LSST 
constrain ildo and wq to 5% and 10%, respectively, 
whereas the strongest SN constraints are derived from 
JDEM, at 10% and 5% for ildo and wq respectively. 
When we combine the probes of lensing and galaxy to- 
mography, SNe, and CMB, the parameter constraints 
improve by a factor of seven in wq and factor of four 
in fldo compared to the constraints from the strongest 
single probe, here weak lensing from LSST. 

The results of the joint analysis don't change sig- 
nificantly when we relax the assumption of spatial 
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FIG. 10. Error ellipses showing degeneracies between cur- 
vature (fifc), early dark energy density (f2e) and sum of 
neutrino masses (^ m^). The curves have the same mean- 
ing as in Fig. |9l Note the strong degeneracy between cur- 
vature and sum of neutrino masses for C MB d ata (black, 
dot-dashed) that was pointed out in Ref. [l52| . This de- 
generacy is broken when information on curvature from 
measures of cosmological distances are included, as shown 
by the solid (red) contour. The constraints in the plane of 
Sle and shown in Fig. [O] 



fiatness. The exception to this statement is Vldo for 
JDEM, which degrades by about a factor of 2 for the 
case where EDE density is fixed (see Table W\ . This 
is because for JDEM Vlda is most strongly constrained 
by SNe measurements, which require a tight bound on 
the curvature. In the joint analysis, the curvature den- 
sity is constrained to 6 x 10^"* of the critical density, 
which is an order of magnitude stronger than solely 
with the CMB temperature and lensing. The ability 
to measure curvature down to this level is an ex citin g 
possibility that has been highlighted previously jl56j |. 
Our constraints on the curvature in the joint analysis, 
and from combining the CMB exclusively with galax- 
ies or weak lensi ng in Fig. I Villi are consistent with the 
results in Refs. |l56l ll57j . even with the introduction 
of early dark energy. 

The sum of neutrino masses is most strongly con- 
strained by the CMB temperature and lensing spectra 
(including their cross-correlation), at the level of 0.20 
eV when the curvature is fixed (Table |IV|. When cur- 
vature is allowed to vary, the sum of neutrino mass 
constraint degrades to 0.22 eV. There is a strong de- 
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TABLE IV. la uncertainties on cosmological parameters from a combination of probes. P denotes CMB T, E, TE modes 
for a Planck survey. K denotes the CMB lensing potential power spectrum and the correlation with the temperature field 
for Planck. L denotes an LSST type survey, whereas J denotes a JDEM type survey, and the superscripts k, g, s, refer 
to weak lensing tomography, galaxy tomography, and supernova measurements, respectively. When we combine more 
than one probe, all relevant cross-correlations between the selected probes are included. Thus, for the case of PKL'^L^, 
all cross-correlations between [T, E, Kc, k, g] are included (see Egn. I48(l . The subscripts with £ < 1000 refer to cutoffs 
of the respective auto-correlations (and all related cross-correlations) ai £ = 1000. At redshifts z = [0, 1,2,3] the early 
dark energy constitutes [0,2.1,8.0, 17.7]% of the overall amount of dark energy (quantified as Qd{z) — i}ui{z) /Qd{z) with 
w = —1 and fie = 0.01). For the case where ^ < 30 modes in CMB polarization data are excluded, we add a prior of 
0.01 on the optical depth. 
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TABLE V. Same as Table [IV] (Planck CMB), for an EDE fiducial cosmology with f7e = 0.01 kept fixed. 



generacy between sum of ne utrin o masses and curva- 
ture as pointed out in Ref. |l52l l. and this is evident 
in the bottom panel of Fig. 1101 However, in the joint 
analysis the eonstraint on curvature (of 6 x lO"** of 
the present critical density) is strong enough to break 



this degeneracy, which will allow for an exquisite mea- 
surement of the sum of neutrino masses. Our fore- 
casts including CMB lensing, weak lensing tomogra- 
phy, galaxy tomography, and supernovae show that 
constraints on the sum of neutrino masses at the level 
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TABLE VI. Same as Table [IV] (Planck CMB), allowing for Qk to vary. We include a CMB prior on the curvature for 
the SN surveys. For the galaxy measurements alone, the bias parameters defined in Sec. Ill Dl as b{z) = bo + biz are 
constrained to a{bo,bi} = {0.16,0.19} for LSST and a{bo,bi} = {0.34,0.43} for JDEM. Our 10% prior on these bias 
parameters (i.e. a{bo, bi} = {0.1, 0.08}) therefore improves the constraints on cosmology. When all probes are combined, 
the bias parameters are constrained to a{bo, fei} = {6.7, 3.5} x lO"'' for LSST and a{bo, bi} = {7.4, 4.6} x lO"'' for JDEM. 
An imposed 10% prior is therefore rendered negligible in the multi-probe scenario. 
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TABLE VII. Same as Table ED without SNe. 



of 0.03-0.04 eV is possible (see Tables |TV] and Ell) • At 
this level, a detection is possible if the neutrino mass 
hierarchy is inverted as Fig. 1111 shows. 

When the early dark energy density and curvature 
are both fixed, the constraint on sum of neutrino 
masses is 0.15 eV from Planck alone. This clearly 
shows that the sum of neutrino masses is substantially 
correlated with early dark energy. This is also appar- 
ent in the right column top panel of Fig. [S] where we 
plot the constraints in the plane of the sum of neu- 
trino masses and EDE density. Thus, as expected, 
ignorance of the expansion history at high redshifts 
can significantly hamper the constraints on the sum 
of neutrino masses from CMB lensing and cosmic 
shear. However, the constraints (even in the presence 
of EDE) can be substantially improved by adding cos- 
mic shear tomography data because the degeneracy 
direction in cosmic shear data is substantially differ- 
ent from that in CMB lensing (see Fig. [9]). 

We note here that our estimate of 0.15 eV when 



curvature and EDE density are held fixed is based on 
the assumption that reionization happened sharply, 
and therefore all of the effects on the CMB can be 
encapsulated in the optical depth parameter. We also 
allow for more complicated reionization histories by 
imposing a floor of = 30 in the polarization data 
and a prior on r of 0.01 [lil. With this setup, we find 
that neutrino mass constraints degrade by about 10% 
(see P£>3oK + cr(r) cases in Tables irvlVllV| . 

We also considered a case assuming sharp reioniza- 
tion where wq is held fixed as is EDE density and 
curvature. In this rather optimistic setup, we obtain 
an uncertainty of about 0.14 eV for the sum of neu- 
trino masses with Planck lensing. We note that this 
is about 20% less stringent than the number quoted 
in Ref. [1^ , some of which can be traced to their dif- 
ferent fiducial cosmology and larger sky fraction (75% 
vs 65%). 

We checked that the cross-correlation of Kc and T 
(i.e. without KcKc) doesn't particularly improve the 
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TABLE VIII. Same as Table IVll but only for combinations of CMB with weak lensing power spectrum measurements, 
and CMB with galaxy power spectrum measurements. 
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0.0038 


0.066 


0.0069 


0.0031 


0.24 


0.0019 


0.23 


0.12 


0.0053 


0.0020 




0.017 


0.0037 


0.0037 


0.059 


0.0066 


0.0032 


0.23 


0.0017 


0.20 


0.13 


0.0052 


0.0015 


^^■'?<1000''^<1000 + ^fe 


0.027 


0.0063 


0.0053 


0.082 


0.0076 


0.0048 


0.24 


0.0022 


0.24 


0.18 


0.0053 


0.0020 


PT^T T ^ T s 4- O. 

-^'^^f<1000^<1000 ~<~ "'■k 


0.012 


0.0029 


0.0033 


0.056 


0.0064 


0.0032 


0.23 


0.0017 


0.20 


0.13 


0.0052 


0.0014 


PK 1" 1® 1^ 4- Oi 


0.011 


0.0030 


0.0047 


0.079 


0.0074 


0.0047 


0.24 


0.0022 


0.23 


0.18 


0.0053 


0.0020 


PKL-^LS + nk-Qe 


0.0094 


0.0023 




0.036 


0.0048 


0.0023 


0.22 


0.0016 


0.19 


0.091 


0.0050 


0.0014 


PRJ-^JS + ^fc - He 


0.015 


0.0037 




0.043 


0.0066 


0.0031 


0.24 


0.0020 


0.23 


0.11 


0.0053 


0.0021 


PKL-^LSL^ + fife - He 


0.0081 


0.0022 




0.036 


0.0048 


0.0023 


0.22 


0.0016 


0.19 


0.091 


0.00.50 


0.0014 


PKJ'^JSJ'' + Hfe - f7e 


0.0081 


0.0026 




0.042 


0.0065 


0.0031 


0.24 


0.0019 


0.22 


0.11 


0.0053 


0.0020 


PKL^^jgggL^^jggg + fife — 


0.012 


0.0034 




0.036 


0.0052 


0.0032 


0.22 


0.0017 


0.18 


0.10 


0.0051 


0.0014 


PKJ^^jggg J^^jggg + fife " fig 


0.021 


0.0059 




0.044 


0.0067 


0.0045 


0.24 


0.0022 


0.22 


0.14 


0.0053 


0.0020 


PKL^^jgQgL^^jgggL^ + fi fe — fi E 


0.0092 


0.0028 




0.036 


0.0052 


0.0032 


0.22 


0.0017 


0.18 


0.10 


0.0051 


0.0014 


PKj^^JQQQ j^^JQQQ J" + fife — fie 


0.0084 


0.0030 




0.043 


0.0067 


0.0044 


0.24 


0.0021 


0.22 


0.14 


0.0053 


0.0020 



TABLE IX. Same as Tables IVll and IVIll without any cross-correlations between the observables (except for tomographic 
ones within a survey). 



constraint on the sum of neutrino masses, but it im- 
proves the constraint on the late time dark energy 
density and EOS by about 50% compared to temper- 
ature and polarization data alone. 

Compared to the CMB with lensing, the LSST 
constraint on the sum of neutrino masses is moder- 
ately weaker for both galaxy tomography (by 50%) 
and weak lensing tomography (by 25%), whereas the 
JDEM constraint on the sum of neutrino masses is a 
factor of three weaker for galaxy tomography and a 
factor of two weaker for weak lensing tomography. In 
analyzing all the probes allowing for EDE density to 
vary, we find a one-sigma constraint of 0.03-0.04 eV for 
LSST and JDEM (see Table |IV|. Tables EIIl and |IX] 
show that the strength of JDEM in this joint analysis 
is owed partly to the superior SN constraints on the 
dark energy EOS (which helps in breaking degenera- 
cies) and partly to the cross-correlations between the 
temperature, lensing, and galaxy fields. 

Given the recent 2a indications of extra relativistic 
species in the universe (e.g. A'ofr = 3.86 ± 0.42 with 



the datasets WMAP7-|-SPT4-HST-fBAO for a mini- 
mal cosm ology jl58j ; for extended parameter spaces 
see |l59{ ). our understanding of N^s depends criti- 
cally on the ability to distinguish its signatures from 
other cosmological parameters. Fortunately, even for 
our extended parameter spaces, Planck alone could 
determine the possible existence of extra relativistic 
species in the universe at the 4(T level, as seen in Ta- 
bles |V] and |Vll The constraints on iVoff are only weakly 
correlated with the constraints on ^ m^, as shown in 
Fig. [121 When Planck is combined with JDEM or 
LSST, this allows for a more than lOa detection of 
additional light degrees of freedom in our universe. 

Turning our attention to systematics. Table IIVI 
shows that the different choices for nonlinear cutoffs 
in ^max S [1000, 3000] do not significantly affect the 
combined constraints. We can therefore throw out 
highly nonlinear scales {£ > 1000) without signifi- 
cant loss in the parameter constraints shown in this 
work. This is primarily due to the combination of 
cosmological probes and to a lesser extent, due to the 
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FIG. 11. Bands show the allowed values for the sum of 
neutrino masses as a function of the lightest neutrino mass 
eigenstate. The branch with higher values is for the in- 
verted mass hierarchy. We have used the constraints on the 
neutrino mass squared differences from the Particle Data 
Group flS^ to create bands of the allowed regions. The 3 
horizontal lines show the la Fisher matrix error estimates 
for the cases PK + fi^ (0.22 eV, Table IVl)) . PK - (0.15 
eV, Table El) and PKL^L^L" + S^^ (0.04 eV, Table ED). 



cross-correlations we have introduced. For compari- 
son, considering only one probe, the constraint degra- 
dation could be up to a factor of 3 in dark energy and 
neutrino mass parameters ( Table ITV]). This insensitiv- 
ity to nonlinear scales is maintained even without SN 



measurements (Tabic IVlTl) . 

Comparing Tabic IIVI with Table I VII i.e. compar- 
ing constraints for a flat universe with a universe 
that allows for a possible non-zero curvature, shows 
the CMB temperature, polarization, and lensing con- 
straints on fifc from Planck improve by an order of 
magnitude when accounting for future weak lensing, 
galaxy clustering, and supernova measurements from 
LSST or JDEM. The SN measurements are, as ex- 
pected, most sensitive to the curvature prior, exhibit- 
ing significant degradations in parameters across the 
board. This motivates the consistent use of a CMB 
curvature prior for SN measurements. The parame- 
ter that the SN observations with a JDEM-like sur- 
vey measure well is the dark energy EOS (the relative 
weakness of SNc from LSST is explained in Sec. IIIip . 
Tables EI] and Ivn] show that the neglect of SN obser- 
vations has < 10% impact on the cosmological con- 
straints for LSST, whereas for a JDEM survey the 
dark energy parameters degrade by up to 50% with- 
out SN measurements. The same features are found 
in the scenarios with an fmax = 1000 nonlinear cutoff 
in multipole space. 

Tables IVII and IIXI show that most parameters im- 
prove by 20% — 40% when cross-correlations are in- 
cluded for LSST. The cross-correlations have a greater 
impact for JDEM than for LSST. In particular, the 
early dark energy density, baryon and CDM densities, 
and sum of neutrino mass constraints improve by up 
to factor of 2 with cross-correlations, the curvature 
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FIG. 12. Error ellipses showing degeneracies between the 
sum of neutrino masses C^rriv) and the effective num- 
ber of neutrinos (A'cii) in a flat universe, where we have 
marginalized over the early dark energy density. The 
curves have the same meaning as in Fig. [5] Even for 
our extended parameter spaces, the CMB temperature as 
measured by Planck will be able to determine the possible 
existence of extra relativistic species, while the combina- 
tion of experimental probes is extremely helpful in pinning 
down the sum of neutrino masses. 



density constraint improves by factor of 3 (also for 
LSST), while the improvement in the other parame- 
ters are at the same 20% -40% level as for LSST. Nat- 
urally, when we weaken the parameter constraints by 
imposing £niax = 1000 in galaxy and lensing tomog- 
raphy, the impact of the cross-correlations increases 
somewhat (e.g. an additional 20% in early dark energy 
and sum of neutrino masses for JDEM) . These quanti- 
tative comparisons between a full covariance and one 
without cross-correlations largely hold true indepen- 
dently of the SN sample. In the final analysis with 
all cross-correlations, both JDEM and LSST are ex- 
pected to provide similar parameter constraints even 
though the individual constraints and systematics are 
different. 

We now turn attention to the impact of keeping 
EDE fixed. In Table|V]we explore the case where EDE 
exists but is kept fixed for the Fisher matrix analysis. 
We expect this to mimic the analysis where the exis- 
tence of EDE is neglected. We explore biases resulting 
from such a scenario in the ensuing section. We expect 
the constraints in this limit (of keeping fig fixed) to 
match the constraints in a non-EDE cosmology, i.e., 
the standard case with ftg = 0. Table |V] shows that 
the CMB temperature and polarization constraints on 
cosmology arc nearly unaffected (< 10% differences) 
by the removal of EDE. The removal of EDE improves 
the sum of neutrino masses by 30%, and most other 
parameter constraints by < 10% when CMB lensing 
is included. 

If we turn to the non-CMB probes, we find that sig- 
nificant degeneracies with fig exist for the dark energy 
parameters, sum of neutrino masses, and the param- 
eters determining the normalization, spectral index. 
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Probe 






wo 




ne 




ris 


d\n k 


IQi^Alj 








T 


Uk 


P + Qk 






0.21 


0.042 


0.0041 


0.32 


0.0042 


0.0049 


0.15 


0.0017 


0.082 


0.14 


0.0031 


0.0069 


PK + Qfe 






0.092 


0.024 


0.0038 


0.079 


0.0040 


0.0044 


0.14 


0.0013 


0.077 


0.11 


0.0029 


0.0024 


PKL'^LfL^ 






0.0062 


0.0015 


0.0017 


0.027 


0.0035 


0.0015 


0.13 


0.00074 


0.064 


0.052 


0.0027 


0.00046 


PKJ«J9J = 






0.0065 


0.0019 


0.0020 


0.029 


0.0037 


0.0019 


0.13 


0.00088 


0.066 


0.058 


0.0028 


0.00061 


PKL^^JQQQ 


T S T S 

«<1000 




0.0065 


0.0016 


0.0018 


0.028 


0.0036 


0.0024 


0.13 


0.00076 


0.066 


0.061 


0.0028 


0.00051 


PK T'^ 


tS ts 




0.0067 


0.0020 


0.0021 


0.030 


0.0038 


0.0029 


0.13 


0.00097 


0.067 


0.069 


0.0028 


0.00065 



TABLE X. Same as Table |Vll except for a 2m EPIC survey in lieu of Planck. 



Probe 




WO 


QdO 




ris 


dris 
din k 


1010A2j 








T 


PKL^LS 




0.83 


0.45 — 


0.23 


0.41 


0.091 


0.37 


0.40 


0.061 


0.41 


0.29 


PKJ«J9 


- He 


2.4 


1.5 — 


0.17 


0.35 


0.27 


0.73 


0.34 


0.029 


0.39 


0.48 



TABLE XL We present the fractional bias \S{pi)\/a{pi) (see Ean. l49|l on cosmological parameters for the full combination 
of lensing, galaxy, and CMB surveys up to ^ = 2000. The letters [L, J] in the first column denote [LSST, JDEM] (precise 
explanation of notation in caption of Table HV)) . The parameter biases are comparable for Planck and EPIC. In this 
table, we assume a flat fiducial universe with f2e = 0.01, and determine the offset in the parameter estimates for the case 
where EDE is not accounted for in the analysis. Naturally, the bias is larger for a universe with a larger fraction of EDE. 



and its running of the matter power spectrum, as ex- 
pected from Fig. [2] For instance, for weak lensing with 
both LSST and JDEM, the power spectrum normal- 
ization constraint improves by more than a factor of 
two without EDE, the present DE EOS improves by 
70%, the rmming by 50%, the sum of neutrino masses 
and spectral index by 30%, and the present DE den- 
sity by 20%. When wc throw out nonlinear scales 
(^max = 1000), the constraint on the sum of neutrino 
masses only improves by about 5%. This is because 
the sum of neutrino masses derives its degeneracy with 
the EDE density from small angular scales, as also il- 
lustrated by the respective parameter derivatives in 

Fig.m 



To end this section, we consider the constraints 
achievable with a future C MB missio n such as the 2m 
EPIC (EPIC-2m) proposal [loS[lOl| in Table B The 
properties of this probe are listed in Table Hill With up 
to 40% increased resolution and factor of thirty lower 
noise in an individual band than Planck, the EPIC-2m 
survey in Table |X] shows improved constraints across 
the board compared to Planck. The joint observa- 
tional constraints, with EPIC-2m in lieu of Planck, 
also show significant gains - about 30% improvement 
in the DE constraints, about 40% improvement in the 
sum of neutrino masses constraint, and up to a factor 
of two improvement in the other parameters. 



In SN measurements, we find the strongest degener- 
acy between the present dark energy density and the 
early dark energy density, which degrades the con- 
straint on the present DE density by factor of five in 
LSST and factor of seven in JDEM. In the joint ob- 
servational analysis of prospective cosmological con- 
straints, we also find that our ignorance of EDE sig- 
nificantly degrades our constraints on the other cos- 
mological parameters. As an example, when keep- 
ing fie fixed in Table |Vl the DE EOS constraint of 
cr(wo) ~ 3 X 10-3 (for both LSST and JDEM) has 
improved by a factor of 3 for LSST and factor of 2 for 
JDEM. In Table we find these degradations would 
have been stronger if SN observations were not con- 
sidered in the analysis. We further note that this 0.3% 
constraint on the EOS may be translated into roughly 
10 mo del- independent redshift bins each constrained 
at the 1% level (since = ^-^(u;) (ill, 

where Wi are binned EOS parameters and w refers to 
a constant EOS). Thus, if the dark energy EOS varies 
at late times, there is significant promise in detecting 
its variation with our considered surveys. 



C. Dark Energy Bias 

We have analyzed the expected constraints on early 
dark energy from use of weak lensing tomography, 
galaxy tomography, SNe, and CMB. If we live in a 
universe with EDE but do not account for this in the 
analysis of the data sets of these future probes, we will 
compute cosmological parameter constraints that are 
not only overly optimistic, but the best estimates of 
the parameters will also be shifted. 

We estimate here the parameter bias that would 
arise from assuming a fiducial universe with fig = 
0.01, but where EDE is not accounted for in the con- 
straint analysis^ The bias in each parameter is given 

by [Mlieffln' 
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M X Tr 



P^-l (^^l ^-l 



(49) 



where 5Ci accounts for the difference between the con- 
vergence spectra in a universe with EDE and a uni- 
verse without EDE. 
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Table IXII presents the cosmological parameter bi- 
ases for the LSST and JDEM surveys combined 
with Planck (comparable biases when combined with 
EPIC). For LSST, we find less than 50% bias in most 
parameters except for wq (about 80%). Biases when 
including JDEM instead of LSST are close to a fac- 
tor of 2 times the expected constraints for the present 
DE density and EOS. We note here that these esti- 
mates for the biases assume a fiducial cosmology with 
Oe = 0.01. The biases would be larger in a universe 
with a larger fraction of EDE. If we are to avoid sig- 
nificant bias in future joint analyses of cosmological 
parameters, we must account for the uncertain high 
redshift nature of the dark energy component. 



IV. CONCLUSIONS 

It is conceivable that in order to resolve the na- 
ture of the force behind the late-time acceleration 
of the universe, we need a deeper understanding of 
how the universal expansion and growth of struc- 
ture was affected by the dark energy at earlier times. 
We have chosen an approach that combines probes 
with varying degrees of sensitivity in different parts 
of cosmological parameter space, and we have in- 
cluded the cross-correlations between different probes. 
We have performed a joint Fisher matrix analysis of 
prospective measurements of weak lensing tomogra- 
phy, galaxy tomography, SNe, and Planck CMB. 

Our results show that the best possible constraints 
on [wo, ildo, fie, Y^m^{eV), n. 



dris 



inlOA2 O 7,2 



lO^ilbh'^, Nee, T, rik] will be at the level of [0.0085, 
0.0019, 0.0023, 0.038, 0.0051, 0.0021, 0.20, 0.0013, 
0.13, 0.090, 0.0044, 0.00056] for a wide survey like 
LSST, and [0.0076, 0.0021, 0.0025, 0.036, 0.0055, 
0.0024, 0.19, 0.0012, 0.14, 0.090, 0.0045, 0.00075] for a 
deep survey like JDEM. These constraint assume that 
photometric redshift and other systematic uncertain- 
ties in weak lensing can be controlled to better than 
0.1%. Naturally, these constraints improve in a uni- 
verse without EDE. For constraints on our extended 
paramet er sp ace with present cosmological datasets, 
see Ref. [l59|. 

Our main findings are summarized below. 

• The dominant constraint on the fraction of dark 
energy at early times comes from the CMB (in 
particular TT and to some extent Tkc). Con- 
straints from Planck are expected to be within 
a percent of the critical density. The next best 
probe of EDE is weak lensing tomography, fol- 
lowed by galaxy tomography, and then SNe. 
When all 6 observables (CMB temperature T, 
CMB polarization E, CMB lensing Kc, weak 
lensing shear k, galaxy counts g, SNe s) and 
all relevant cross-correlations (see Eqns. ITTIHS]) 
are analyzed in a combined setting, the con- 
straint on fie improves by a factor of four over 



the CMB constraint. Interestingly, due to the 
cross-correlations, the combined constraints are 
equally strong for LSST and JDEM (k, g, s). 

• When analyzed together with early dark energy, 
we find that the sum of active neutrino masses is 
most strongly constrained by the CMB lensing 
potential power spectrum, at the level of 0.22 
eV for Planck (0.20 eV when flatness assumed) 
and 0.08 eV for Epic (0.07 eV when flatness as- 
sumed). A combined analysis of Planck, and 
LSST or JDEM, shows that future constraints at 
the level of 0.04 eV are possible. The CMB lens- 
ing constraints improve by 30%, and the com- 
bined constraints by 60%, when EDE is not al- 
lowed to vary. Our results suggest that these 
constraints are not significantly affected by our 
ignorance of the reionization history, but more 
detailed work on this issue is necessary. 

• The additional degree of freedom from the early 
dark energy density degrades our ability to mea- 
sure late-time dark energy. We find that the 
present DE density can be measured to 0.2% of 
the critical density and and equation of state to 
about 0.01, which is a factor of roughly 2 and 
3 degradation, respectively, in the constraints 
compared to the case when EDE is fixed. 

• Our analysis suggests that throwing out nonlin- 
ear scales {i > 1000) may not result in signifi- 
cant degradation in future parameter measure- 
ments when multiple cosmological probes are 
combined. Including cross-correlations improves 
parameter constraints on dark energy density 
and sum of neutrino masses by up to a factor of 
2 when these nonlinear scales are not included. 

• The curvature of the universe can be constrained 
to 6 X 10~^ of the critical density from CMB 
temperature and lensing alone, and improved by 
an order of magnitude in the joint ana lysis in 
agreement with the results of Ref. 1571 ]. Mea- 
surement of the curvature of the universe with 
Planck will be good enough that weak lensing 
tomography, galaxy tomography, and supernova 
measurements will not be limited by our igno- 
rance of the curvature. 

• Even a modest 1% of the critical density at high 
redshift in dark energy, if not accounted for, 
shifts cosmological parameters by 1-2 a. There- 
fore, it is crucial for measurements of the under- 
lying cosmology that we avoid prejudices about 
the energy content of the high redshift universe. 

We have shown that degeneracies between cosmo- 
logical parameters, in particular between early dark 
energy, curvature and sum of neutrino masses, can be 
effectively broken by a joint analysis of weak lensing 
tomography, galaxy tomography, SNe, and the CMB 
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(temperature, polarization, lensing). Our analysis in- 
cluded the cross-correlations between these different 
probes. We find that simultaneous measurements of 
dark energy density at present and at high redshift 
with a precision of 0.2% of the respective critical den- 
sities, present dark energy equation of state with a 
precision of 0.01, curvature with a precision of 0.06% 
of the present critical density, and sum of neutrino 
masses with a precision of 0.04 eV are possible. 
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